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We report here on the application of the perturbative renormalization-group to 
the Coulomb gauge in QCD. We use it to determine the high-momentum asymptotic 
form of the instantaneous color-Coulomb potential V(k) and of the vacuum polariza- 
tion P(k,k4). These quantities are renormalization-group invariants, in the sense that 
they are independent of the renormalization scheme. A scheme-independent defini- 
tion of the running coupling constant is provided by k 2 V(k) = xog 2 (k/ A cou i), and of 
a s = 9 ( fc /^ cou ') ; where xq = 11 ^j 2 ^ Nf , and A cou i is a finite QCD mass scale. We 
also show how to calculate the coefficients in the expansion of the invariant /9-function 
(3(g) = |fc|-J§: = — (bog 3 + big 5 -\-b2g 7 + ■ ■ .), where al I coefficients are scheme-independent. 
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1. Introduction 

In QCD the Wilson loop is the basic gauge-invariant observable. A rectangular 
Wilson loop of dimension R x T has, asymptotically at large T, the form W(R, T) ~ 
exp[— TVw(R)]- The Wilson potential Vw(R) is the total energy of the state of lowest 
energy that contains a pair of infinitely massive external quarks at separation R. If dy- 
namical quarks are present in the theory, the pair of external quarks polarizes the vacuum 
and extracts a pair of dynamical quarks from it, so that, for R not too small, a pair of 
mesons is formed at separation R, each meson being formed of one external quark and 
one dynamical quark, plus doses of other constituents. In this case Vw(R) represents the 
total energy of the pair of mesons at separation R. Clearly the Wilson potential Vw(R) is 
not a color-confining potential, but rather a residual potential that remains after color has 
been saturated by vacuum polarization. In this respect Vw(R) should be regarded as a 
QCD analog of the van der Waals interatomic potential. Clearly it is not the van der Walls 
force that holds the atom together, but rather the electrostatic Coulomb potential. In the 
present article we shall be concerned with the QCD analog of the electrostatic Coulomb 
potential, which we regard as responsible for confinement of color charge. 

This quantity is the the color-Coulomb potential V(R). Just as the electrostatic 
Coulomb potential is the instantaneous part of the 44-component of the photon propagator 
in the Coulomb gauge, likewise V(R) is the instantaneous part of the 44-component of the 
gluon propagator in the (minimal) Coulomb gauge, defined in Eq. (1.3) below. It is not 
gauge invariant like the Wilson potential Vw(R), but it is a more elementary quantity, 
in terms of which one may hope to understand the dynamics of the underlying theory. 
Indeed we expect that its linear rise (or not) at large R may serve as an order-paramenter 
for the confinement of color-charge, even when V\y {R) is not linearly rising due to vacuum 
polarization. In the present article, we report on a renormalization-group calculation of 
V(R) in the high-momentum regime, more precisely, a calculation of its Fourier transform 
V(|fc|) at large \k\. 

To be sure, confinement manifests itself rather at low momentum. For information in 
this region we have turned to numerical study, and in an accompanying article on SU(2) 
lattice gauge theory without quarks [1], we present a numerical determination of 
and also of the equal-time would-be physical 3-dimensionally transverse gluon propagator 
Dfj(k). The reader may find further references in [1], and a confrontation of the confine- 
ment scenario in the Coulomb gauge with the numerical data. This confinement scenario 
is also discussed in [2], and in Sect. 3 of the present article. 
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The Coulomb gauge is a "physical gauge" in the sense that the constraints are solved 
exactly, including in particular Gauss's law, DiEi = p qu , that is essential for color confine- 
ment. In a confining theory, all physical states are bound states, and the Coulomb gauge is 
the preferred gauge for calculations of bound-states. Binding is provided by the instanta- 
neous color-Coulomb potential V(R) that is treated non-perturbatively, while everything 
else is regarded as a perturbation [3] , [4] . However the Coulomb gauge offers no particular 
advantage for purely perturbative calculations. 

In the present article we shall apply the perturbative renormalization-group in the 
Coulomb gauge. This presents new challenges. However they are amply rewarded by 
particularly strong results. These are a consequence of the fact that in the Coulomb gauge 
the 4-component of the vector potential is invariant under renormalization [2] , 

g A 0A (x) = g r A rA (x), (1.1) 

where the subscripts and r refer to unrenormalized and renormalized quantities. This is 
not true in a Lorentz-covariant gauge. 

The Coulomb gauge is traditionally defined by the condition diAi = 0. This is an 
incomplete gauge fixing because gauge transformations g(t) that depend on t are not 
fixed, and consequently in higher-order calculations one encounters singular expressions 
whose evaluation is ambiguous [5], [6]. No doubt this has been an obstacle to the use 
of the Coulomb gauge in QCD. These difficulties are overcome however by defining the 
Coulomb gauge as the limit of the interpolating gauge characterized by the condition 
diAi + \d4A4 = 0, where A is a real positive parameter [7]. Calculations are done at finite 
A, and at the end one takes the limit A — > 0. The interpolating gauge is a renormalizable 
gauge that may be treated by standard BRST methods. The gauge-fixing term breaks 
Lorentz invariance, but it is BRST-exact. An extension of the BRST operator to in- 
clude infinitesimal Lorentz transformations assures Lorentz invariance for BRST-invariant 
observables [7]. This definition of the Coulomb gauge is sufficient for the purposes of 
perturbation theory, and for the perturbative renormalization-group that is the subject of 
the present article. In fact the only calculations that we will do explicitly (in Appendix 
B) are one-loop, and for these one may set A = directly. However in our discussion of 
the renormalization group we rely on the existence of the limit A — > in every order of 
perturbation theory. 
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[At the non-perturbative level with, say, lattice regularization, the interpolating gauge 
condition is still subject to the Gribov ambiguity. This is resolved by the further specifica- 
tion of the minimal Coulomb gauge, whose lattice definition is given in [1] . Its continuum 
analog would be to first minimize Fhor,A(g) = / d 3 x 9 A 2 i , for each t, with respect to all 
local gauge transformations g(x, t), where 9 A fl = g~ 1 A^ l g + g~ 1 d^g is the gauge transform 
of A. At the minimum, the Coulomb gauge condition diAi = is satisfied for each t, and 
the Faddeev-Popov operator M(A) = —Di(A)di is positive. Next F ve r,A(g) = / d 4 x 9 A A 
is minimized with respect to t-dependent gauge transformations g(t). At the minimum 
/ d xA%{x, t) = C, where C is a constant that satisfies C < 2tvV/T. Here V is the spatial 
3-volume and T is the temporal extent of the Euclidean 4-volume.] 

An immediate consequence of (1.1) is that the 44-component of the gluon propagator 
is also invariant under renormalization. This gives the relations 

D 44 (k, fe 4 , Aqcd) = D QM {k, k A ,g , A) = D rM {k,k±,g r , /i). (1.2) 

(We use the same symbol to designate a quantity and its Fourier transform, but dis- 
tinguish them by their argument x or k.) Here Aqcd is a finite QCD scale, A is the 
ultraviolet cut-off, and \i is the renormalization mass. We have written the finite quantity 
D^k, &4, Aqcd) without a subscript or r because it is a renormalization-group invariant 
in the sense that it is independent of A and \i and of the regularization and renormaliza- 
tion schemes. The finite scale Aqcd, for example Aqcd = Ai att or Aqcd = A™ s , appears 
in D , 44 (fc, &4, go, A) only through its dependence on g = g (A/ Aqcd), and similarly for 
D r ,44(k, k^g r ,[i). The scale is defined within the particular regularization scheme that is 
used for calculations, for example lattice or dimensional regularization. To relate two dif- 
ferent scales, one calculates the same renormalization-group invariant in the two schemes. 

In the Coulomb gauge D±±(x,t) has the decomposition into an instantaneous part, 
proportional to 8{t) and a non-instantaneous part P{x, t) that is less singular at t = 0, 

D 4 4x,t) = V{x)8{t) + P{x,t), (1.3) 

where V(x) is what we call the instantaneous color-Coulomb potential, and P(x, t) is a 
vacuum polarization term. (See Eqs. (2.13)-(2.15) below.) We shall see that V(x) is 
anti-screening, whereas P(x, t) is screening. In momentum space this decomposition reads 

D M (k, k 4 ) = V(k) + P(k,k 4 ). (1.4) 
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The instantaneous part is independent of k 4 . The ^-independent part is defined, without 
reference to a diagrammatic expansion, by 

V(k)= lim D 44 (k,k 4 ). (1.5) 

Here it is assumed that the limit (1.5) exists and is independent of the sign of k 4 . This 
implies that 

P(k,k 4 ) = D 44 (k,k 4 ) -V(k). (1.6) 

vanishes at large k 4 , 

lim P(k,k 4 ) = 0. (1.7) 

The only tool available to calculate these quantities analytically is perturbation the- 
ory. Perturbatively, in d < 4 dimensions, the separation into instantaneous and vacuum 
polarization parts may be made diagram by diagram. Each diagram is either independent 
of k 4 or satisfies (1.7). (See the decomposition (2.13) - (2.15) below.) However in d = 4 
dimensions, individual vacuum polarization diagrams contain powers of ln(Ak 4 ) associated 
with divergences, and the limit (1.5) does not exist in any finite order of perturbation the- 
ory. As a result there is an ambiguity in the separation of the instantaneous and vacuum 
polarization parts in any finite order of perturbation theory. (This may be regarded as a 
mixing of V and P due to renormalization.) However because D 44 is a renormalization- 
group invariant, it follows from the above definition that V(k) and P(k, k 4 ) are separately 
renormalization-group invariants. This will allow us to resolve the ambiguity by use of the 
perturbative renormalization-group. We shall then find the asymptotic form of V(k) at 
large \k\ and of P(k, k 4 ) at large k 4 . 

As a result of a renormalization-group calculation, we find that in an SU(N) gauge 
theory with Nf quarks in the fundamental representation, the instantaneous color-Coulomb 
potential is given by 

k 2 V(k)= x g 2 (\k\/A coul ), (1.8) 

where 

12N , 
X ° = llN-2N f > (L9) 

and A cou i is a finite QCD mass scale specified below. The running coupling constant 
g(\k\/ A cou \) is found as the solution of the renormalization-group equation 

\k\-^r=fi(g), (1.10) 
o\k\ 
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where the /3-function has the weak-coupling expansion 



P{g) = -(bog 3 + hg 5 + b 2 g 7 + ...), (1.11) 

where all coefficients b n are scheme-independent, not just bo and b\. The coefficients bo 
and b± are also gauge-invariant and have their standard values. We shall show how to 
calculate the remaining coefficients perturbatively in the Coulomb gauge. This allows 
a calculation of V(k) to arbitrary accuracy in the high- momentum region. Its leading 
asymptotic behavior is given by 

Pv ^ ] ~ 9 , w^i/a (i-i2) 

26 ln(|/c|/A CO ui) V ' 

We suppose that asymptotic freedom holds, which requires that bo, given explicitly below, 
be positive, bo > 0. This is equivalent to lliV — 2Nf > 0. We observe that xq > 0, and in 
fact xq > 1. 

We have also found the leading asymptotic behavior of the vacuum polarization term 
P(k,k 4 ) at large k± namely, 

pp-(fl.*..»,A)~ 2Mn( r 4 /A- ou , r ^ 

where A^. oul is another finite QCD mass scale, and y = 1 — x = — iiN-2N f < ^- 
negative sign of yo shows that the vacuum polarization term is indeed a screening term. 



2. Formula for D 44 in the Coulomb gauge 

In the minimal Coulomb gauge, the Euclidean partition function is given by the fa- 
miliar Faddeev-Popov formula 

Z(J, J 4 ) = J d 4 A 6(diAi) det(-A^) exp J d 4 x [-(1/4)^ - igoJ^)}, (2.1) 

where is the Yang-Mills field. We have introduced sources J M , in terms of which the 
gluon propagator is given by 

iW* - V»* - sKA^AUy)) = -Z-^^U,, (2.2) 

(We shall frequently suppress color indices.) The subscript G on the integral refers to the 
fact that the integral over A\ r is restricted to within the Gribov region. This has very 



important dynamical consequences at long range that have been proposed as a mechanism 
for the confinement of color charge [8] and [9], and which are substantiated by numerical 
studies of the Coulomb gauge [1]. However in the present article we are interested in the 
asymptotically high-momentum region where asymptotic freedom reigns, and we need not 
specify the Gribov region G. 

The partition function may also be expressed in terms of the first-order or phase-space 
functional integral, by introducing a Gaussian integral over an independent color-electric 
field Ef. This is done in Appendix A, with the result that, in the minimal Coulomb gauge, 
the partition function may be expressed as an integral over the physical, canonical degrees 
of freedom A\ r and E\* with the canonical action 



Z(J, J 4 ) = f dA tr dE tr 
Jg 



(2-3) 

x exp / d 4 x {iEfAf -H- igoJiAf). 



Here 



H = (l/2)(Ef + Bf) (2.4) 

is the classical Hamiltonian density, Bi = <9 2 A 3 — <9 3 A 2 + goA 2 x A 3 , etc., where Af = Af' a , 
(X x Y) a = f^X h Y c , and the color-elect ric field Ei is expressed in terms of the canonical 
variables by solving Gauss's law, 

DiEi = g J 4 , (2.5) 

where Di = di + goAiX is the gauge-covariant derivative. 

Gauss's law is solved by separating the transverse and longitudinal components of Ei 
according to 

Ei = Ef - d i( j>. (2.6) 

Here cf> is the color-Coulomb potential in terms of which Gauss's law reads —Didi4> + 
g Af x Ef = g J 4 , or 

M(A tr )<f) = p coul + g J 4 , (2.7) 

where M(A tr ) = — Di(A tr )di is the 3-dimensional Faddeev-Popov operator, and p cou i = 
— goA\ T x E\ Y is the color-charge density of the dynamical degrees of freedom. If we had in- 
cluded dynamical quarks then, in addition to a Dirac term q(-fiDi+m)q in the Hamiltonian 
density, there would also be a quark contribution to the color-charge density 

Pcoui = Plx + <u = ~9of abc Af ' b Ef > c + g Q Wot a q. (2.8) 
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The solution of Gauss's law is given by 

(j)(x,x 4 ) = J d 3 y M _1 [x,i/,A tr (x4)] (Pcoui + So^X^O- (2.9) 

Because M(A tr ) involves only spatial derivatives, the inverse operator M~ 1 (A tr ) acts in- 
stantaneously in time. The Hamiltonian is given by 

H = J d 3 xH = (1/2) J d 3 x (Ef + B 2 ) 

= (1/2) fd 3 x [ET 2 + W) 2 + Bf] 

J (2.10) 
= (1/2) J d 3 x (ET 2 + B 2 ) 

+ (1/2) J d 3 x d 3 y (pcoui+9oJ4)(x) V(x,y;A tr ) (p coul + g J^)(y). 

Here 

V(x,y;A^) = [M{A^)-\-d 2 )M{A^)-%,y (2.11) 

is a color-Coulomb potential-energy functional, that depends on A tr . It acts instanta- 
neously, and couples universally to color charge including, in the present case, the source J4. 
Its physical significance will be discussed in the next section. 

We now come to the important point. We use Z from (2.3) and H from (2.10), and 
obtain 

^Jlfix) = - 9o J d3z{ Va W;^ tr ) [Pcoui(^4) +^4(^4)] ), (2.12) 
To calculate Do, 44 (^ — y) we apply 5J t^ and we obtain two terms, 

D 0M (x -y) = V (x -y) + P (x - y). (2.13) 

The first, 

V (x - y)5 ab = g 2 ( V ab [x, y; A tr (x 4 )] ) 5(x 4 - y 4 ), (2.14) 

comes from the J2(z, X4) that appears explicitly in (2.12). It is manifestly instantaneous. 
The second term, 

P (x - y)d ab = -g 2 ( (Vp coul ) a (x) (Vp coul ) b (y) ), (2.15) 

where (Vp) a (x) = (Vp) a (x, X4) = j d 3 z V ab [x, z; A tr (x4)]p b (z, £4), represents polarization 
of the vacuum. The expansion of Vo and P up to one- loop is given in Appendix B. The 
significance of these two terms for confinement is discussed in the next section. 



3. Physical Interpretation 

We comment briefly on the physical meaning of the decomposition (2.13) to (2.15). We 
showed in [1] that the Faddeev- Popov operator M(A tT ) is a positive operator in the minimal 
Coulomb gauge. It follows that V(A tr ) = M~ 1 (—d 2 )M~ 1 is also a positive operator. The 
Pq term represents ordinary vacuum polarization that also occurs in QED. The minus sign 
that appears in front of Po indicates that it corresponds to screening. In QED, the color- 
Coulomb potential energy functional V(x, y, A tr ) would be replaced by the electrostatic 
Coulomb potential (4n\x — y|) _1 . In [1] we argued that V(A tr ) is long-range for a typical 
configuration A tr , and presented numerical evidence that this is true. 

According to the confinement scenario in the minimal Coulomb gauge, V(A tr ) is pre- 
dominantly long range because M _1 (A tr ) is long range. This happens because the Gribov 
region is bounded by M(A tr ) > 0, and entropy favors dense population close to the bound- 
ary. The boundary occurs where M(A tT ) has a zero eigenvalue, so it has a small eigenvalue 
close to the boundary. As a result M~ 1 (A tr ) is enhanced close to the boundary where the 
population is dense. The same color-Coulomb potential energy functional V(A tT ) appears 
in both the instantaneous term Vo and the vacuum polarization term P , so both are long 
range. 

The instantaneous term Vo is responsible for confinement, whereas the vacuum polar- 
ization term P causes screening. In any physical process both Vo and P contribute of 
course, and one or the other may dominate. Consider the Wilson loop which is a model of a 
pair of infinitely heavy external quarks. In a theory of pure glue without dynamical quarks, 
it is believed that at long range there is a linearly rising potential Vw(R) in the Wilson 
loop, characterized by a string tension. In this case, according to the confinement scenario 
in the Coulomb gauge, the instantaneous term Vq dominates. However if dynamical quarks 
are present then it is believed that the string "breaks" at some distance because it is en- 
ergetically favorable to produce pairs of dynamical quarks from the vacuum. In this case 
the vacuum polarization term Pq dominates. However in both cases color-charge itself is 
confined. According to the confinement scenario in the Coulomb gauge, that is because the 
instantaneous term Vq(R) is long range and presumably linearly rising even when Vw(R) 
is not. Indeed it is precisely the long range of Vq(R) that makes it energetically favorable 
to produce dynamical quark pairs from the vacuum, thereby causing color confinement. 
Thus in the Coulomb gauge the instantaneous part provides an order-parameter for the 
confinement of color charge even when dynamical quarks are present. 
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As noted, the term Vo is instantaneous. It is easy to see that in each order of pertur- 
bation theory Vq{x — y) is the sum of diagrams in which the points x and y are continuously 
connected by instantaneous free gluon propagators, as illustrated in Fig. la. We now show 
that the diagrams contributing to Pq{x — y) have no instantaneous parts. Indeed the free 
propagators of A tr and E tr are given by 

(AfAf) = P ij (k)/(P + kl) 

{E?Ef) Q = (k%j - hk^/ik 2 + k\) (3.1) 
(E?Af) = P tJ Ck)h/(k 2 + kl), 

where Pij(k) = (5ij — kikj) is the 3-dimensionally transverse projector. These propagators 
vanish in the limit k 4 — > oo, as do the quark free propagators. From the structure of 
the vacuum polarization term (2.15), as illustrated in Fig. lb, one sees that the diagrams 
contributing to it must somewhere be connected by the dynamical propagators (3.1), so 
Po has no instantaneous part. 2 This is apparent in the one-loop calculation presented in 
Appendix B. 

4. A technical difficulty 

In Eq. (1.5) we have given a definition of the color-Coulomb potential which is inde- 
pendent of perturbation theory. However since its defining property is false in any finite 
order or perturbation theory, we must address the problem of how to calculate it. 

The perturbative expansion of -Do, 44 = Vo + Po begins with 

V (\kl g , A) = g 2 (k 2 )' 1 + O(g 4 ) 
Po(\klk 4 ,g ,A) = 0(g^. 

2 It is essential to use the phase-space integral with first-order action to separate out the 
instantaneous diagrams in this way. For one could also integrate out A4 in the second-order 
Lagrangian formalism. However if one did that, then Af would appear at vertices instead of the 
independent field variable Ef . The free propagator (A* r A* r )o = Pij(k)k\l (k 2 + k\) does not vanish 
in the limit — > 00, but rather has the limit (A* r A* r )o — > Pij(k) which has an instantaneous 
non-local part. 
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The Feynman integrals for these quantities are derived to one loop order in Appendix B. 
The integral for Po(|&|> &4> <7o> A) is relatively complicated, however it simplifies consider- 

— * 

ably for fc 4 » \k\, as discussed in Appendix B, and one obtains to one- loop order for 
SU(N) gauge theory, 

k 2 V (\k\,g , A) = g 2 + g*[ Vll ln(A/|fc|) + i; 10 ] 

(4.2) 

k 2 P (\k\, fc 4 , go, A) = g^\p u ln(A/fc 4 ) + Pio + o(|fc|/fc 4 )], 
where A is an ultra-violet cut-off, lim^o o(x) = 0, and v nm and p nm are constants, with 

8N , s 

^"(Jf^ + N' < 44 > 

Here we have included the vacuum polarization contribution of Nf flavors of quark in the 
fundamental representation which, to this order, is gauge-independent. These coefficients 
are also found in a Hamiltonian formalism [10]. 

We obtain asymptotically for fc 4 » \k\, for D^ 44 = Vq + Pq s , 

k 2 D^ M (\k\,k 4 ,g ,A) = g 2 +g^[v 11 \n(A/\k\) + Pll ln(A//c 4 ) + uio +Pio]- (4-5) 

Observe that pio contributes a 1/k 2 term to -Do, 44- Consequently we cannot assert that 
the instantaneous part of -Do, 44 in d = 4 dimensions, namely the color-Coulomb potential 
V is given by the sum of instantaneous diagrams, as it is in d < 4 dimensions, 

V(\k\,A Q cD)^V (k,g ,A). (4.6) 

We shall see that Vq and Pq separately are not renormalization-group invariants, and this 
will be the key to an unambiguous separation of the instantaneous part V. This implies 
that the coefficients that appear in Vo and Po ma Y be scheme dependent. In particular, 
vio and pio are scheme-dependent. 

Another way to see that the separation cannot be made diagram by diagram is to 
recall that in higher-order, calculations must be done in the interpolating gauge, with 
gauge condition diAi + A<9 4 A 4 = 0. For finite A no diagram is instantaneous, and a 
separation criterion is required for the limiting expression at A = 0. 

As a first step we note that in one-loop order -Dq^ 44 (|A;|, A; 4 , <7o> A) at large fc 4 conve- 

— * 

niently separates into a sum of terms that depend respectively on \k\ and on fc 4 . To get the 
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instantaneous part, namely the part that is independent of k±, we simply delete the term 
that depends on k^, namely pu ln(A//c4). However the separation of the constant term is 
ambiguous, as it is in each order. We conclude that to one-loop order, the /^-independent 

— * 

part of D 0M (\k\, k 4 , g Q , A), is given by 

k 2 V(\k\,g , A) = g 2 x + <7 >n ln(A/|fc|) + Xl ], (4.7) 

where xq and x\ are as yet unknown constants. 

The form of the higher-order contributions to Pq is not known. However we shall make 
the simplest assumption namely that the same procedure may be effected to arbitrary 
order, so that only constant terms in each order in the expansion of k 2 V{\k\, g , A) are 
ambiguous, 

oo n 

£V(|£|,<7o,A) = g 2 x + Y.9o n+2 i*n + £ v nm \n m (A/\k\)}. (4.8) 

n=l m=l 

Here the coefficients v nm are calculated from the perturbative expansion of Vo, but the x n 
are a set of as yet unknown constants. (In case there are additional ambiguous terms in 
higher order in the separation of the vacuum polarization part Po from -Do,44> they may 
also be determined using the renormalization-group, because it also restricts the v nm .) 

From its definition as the large-/c4 limit of D44, it follows that color-Coulomb po- 
tential is a renormalization-group invariant. So when V(|fc|) is expressed either in terms 
of unrenormalized or renormalized quantities, it is independent of A and [i (and of the 
regularization and renormalization scheme), 

V(\k\) = V(\k\,go(A/ Aqco), A) = V(\k\, qM^qcd), /*)• (4-9) 

Consequently we may set A = and fx = \k\, which gives 

V(\k\) = V(\klg (\k\/A QC D), \k\) = V(\klg r (\k\/A QCD ), \k\). (4.10) 

(Once the functional dependence of -Do, 44 or Vo on the cut-off A is determined - it is a 
polynomial in In A in each order of perturbation theory - then A may be assigned any 
finite value.) Thus we may set A = \k\ in (4.8), which gives a simple expansion in terms 
of the unknown coefficients x n , 

00 

k 2 V(\k\) = xo gl{\k\/A QCD ) + Y,*n g 2 n+2 (\k\/A QC D). (4.11) 

n=l 
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5. Renormalization-group in the Coulomb gauge 

The unrenormalized coupling constant go = go(A/ Aqcd) is a function of the cut-off 
A, determined by the flow equation 

Adg /dA = p (g ), (5.1) 

where Aqcd is a finite physical QCD mass scale. The /3-function has the expansion 

&(0o) = -(&o0o + &i0o + &20o + --O • ( 5 - 2 ) 

In general, the coefficients b n are both scheme and gauge dependent, except for 6 and 
b± that are scheme and gauge independent. Similarly the renormalized coupling constant 
9r = 9r(p-l Aqcd) is a function of the renormalization mass determined by 

fi dg r /dii = (3 r (g r ). (5.3) 

In the Coulomb gauge -Do, 44 has no anomalous dimension coming from multiplicative 
renormalization. It therefore obeys the simple Callan-Symanzik equation 

[Ad/dA + (3 {g )d/dg Q ]D 0M {\k\M,g ,A) = 0. (5.4) 

As a result, in the Coulomb gauge, the /3-function may be obtained from the propagator 
Do, 44 - in fact only D^ 44 is needed - whereas in covariant gauges a calculation of the 
vertex function is necessary. Indeed from the last equation we have 

AdD QM /dA 

Po #o = — ^7. jr. — , 5.5 

aD ,44/%) 

where -Do, 44 = -Do,44(fI) ^4><7o> A). This holds identically for all values of |fc|, and A, so 
we may set k± to an asymptotically large value and obtain 

AdDf,JdA 

From the one-loop expression (4.5) for -Dq^ 44 , we obtain the first coefficient of the /3-function 

bo = ^(v 11 +p 11 ). (5.7) 
12 



From (4.3) and (4.4) we obtain the standard expression 

b > = l£r(T N -l N >)< < 5 - 8 > 

without calculating any vertex function. All coefficients b n may be calculated in this way. 
We take the large fc 4 limit of the Callen-Symanzik equation (5.4), and observe that 

— * — * 

V(\k\,go,A) = lim^^oo -Do, 44(1^1, 9o, A) satisfies the same Callan-Symanzik equation, 

[Ad/dA + po(go)d/dg }V(\k\,g , A) = 0. (5.9) 

Since go = go(A/ Aqcd) is a solution of the flow equation (5.1), the Callan-Symanzik 
equation yields 

dV(\klg (A/A QCD ),A) _ Q (510) 
dA 

which is again the statement that V(\k\, go(A/ Aqcd), A) is independent of A. 



6. Renormalization-group to the rescue 

We shall use the renormalization-group and our knowledge of the /3-function to de- 
termine the unknown constants x n . From the Callan-Symanzik equation (5.9) for V we 
have 

avm, !l0 ,A)/am= - Aav f; 9 ' , : A)/aA . 

This is an identity that holds for all A, and we may simplify it by setting A = \k\, 

dV(\k\,g ,A)/dg \ A=lUl = ^ . (6-2) 

To find all the x n we substitute the expansion (4.8) on the left and right hand sides. We 
also expand 

oo 

PoHgo) = ^ " 3 (1 + J2 c p9o P ^ ( 6 - 3 ) 

P =i 

where, from (5.8), we have c\ = —bi/b etc. The derivative with respect to A on the 
right-hand side of (6.2) kills the constants x n , and we have 

oo oo oo 

2(n + l)g 2 n+1 x n = b^g^ 3 ^ c p g 2p £ g 2 m+2 v ml , (6.4) 

n=0 p=0 m=l 
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where cq = 1. Equating like powers of g$ we obtain 

n+l 

x n = [2(n + 1)6 ] _1 ^ Cn— m+l"ml j 

(6.5) 

m=l 

and in particular 

x = (2b )~ 1 v 11 

(6.6) 

Xi = (46 ) (vnci + v 2 i). 

Thus xq is found from vu and &o which require 1-loop calculations, and x\ requires 2- 
loop calculations. Usually the renormalization-group is used to determine higher-order 
logarithms from lower order terms. Here instead we have used it to consistently determine 
an unknown lower-order constant from a known higher-order logarithm. From (4.3) and 
(5.8), we obtain 

12AT / 
X0 = UN-2N/ (6 ' 7) 

It will be convenient in the following to factorize the coefficient xq out of the expansion 
(4.11) for V, and we write 

oo 

k 2 V{\k\) = x [ 9l(\k\/^ Q CD) 9 2 o n+2 (\k\/A QC D) }■ (6.8) 

n=l 

where x' n = x u /xq. Remarkably, the leading term is not simply <7q, as one would expect 
from (4.2), but rather x g 2 = ytw^n] 9o- 

The leading asymptotic form of V may be found from the expression (4.8) to order 
<7q, with vu = 2boXo, namely 

k 2 V(\k\) = x g 2 { 1 + g 2 [2b ln(A/|fc|) + x[}}, (6.9) 

which to this order may be written 

k 2 V(\k\) = x { go 2 ~ [2b ln(A/|£|) + x[] }"\ (6.10) 

where 

^ = - = \{--h- (6-11) 

1 x 2 y v n b J V ' 

From the asymptotic form of g^ 2 ~ 26 h^A/Agci)), this gives 

k 2 V{\k\) ~ x [ 2boH\k\/A QC D)-x' l }- 1 (6.12) 
14 



which is valid for large \k\. Here we have introduced the new mass scale A cou i characteristic 
of the Coulomb gauge. It is related to the scale Aqcd used in the scheme by which the 
/3-function was defined according to 

Acoui = exp (^M Aqcd- (6-14) 

Normally, a ratio such as A cou i/Aqcd can be determined from a 1-loop calculation, but 
here a 2-loop calculation of v 2 i, the coefficient of a 2-loop logarithm that appears in Eq. 
(6.11) would be required. For just as the determination of the constant xq, of 0-loop order, 
requires the 1-loop calculation which we report here, similarly, to determine the constant 
x[, of 1-loop order, would require a 2-loop calculation. In both cases the reason is the 
ambiguity in the separation of the instantaneous part V which, as noted above, cannot 
be determined from individual perturbative diagrams. Instead it requires the perturbative 
renormalization group to determine the n-th order constant from the (n + l)-th order 
logarithm. 

Expression (6.13) for V^(|A;|) exhibits asymptotic freedom and, with x$ positive, indeed 
xq > 1, the instantaneous part of -D44 is anti-screening. 



7. Asymptotic form of the non-instantaneous part 

We may also determine to this order the asymptotic form of the non-instantaneous 
part P = -Do, 44 — V ■ It is also a renormalization-group invariant. If we set to an 
asymptotically high value we obtain 

P as (|£|, fc 4 , go, A) = D? M (\k\, fc 4 , go, A) - V(\k\,g , A). (7.1) 
From (4.5) and (4.7) this gives to one- loop order 

k 2 P as (\k\,k 4 ,g ,A)=g^y +gt\piiH^/k4) + yi}, (7.2) 

where, by (4.4), (5.7) and (6.6), 

yo = l-x = (26o) _1 (26 - vu) = (26 )"Vi 

N + 2N f (7.3) 
llN-2N f ' 
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and yi = pm + vm — x\. This gives 

k 2 P as (\kl h,g , A) = y g 2 {1 + g 2 [2b ln(A/fc 4 ) + y[}, (7.4) 

where y[ = yi/yo- By the reasoning that leads to (6.13), we obtain to this order the 
asymptotic expression 

PP-(|£|,*.,a,,A)~ 2b(M !° i/A , QcD y (7.5) 

where A'q CD is another finite QCD mass scale like A cou i. It would take also take a 2-loop 
calculation to determine Aq CD / Aqcd- Since y < is negative, P is indeed screening. 



8. Invariant color charge 

In this section we show how to calculate V(|/c|) to arbitrary accuracy in the ultra-violet 
region. We define a new running coupling by constant g = g(\k\/ A cou \) by 

k 2 V(\k\)=x g 2 (\k\/A coul ). (8.1) 

Because V(|fc|) is scheme-independent, g(\k\/A cou \) is also. From (6.8) we have 

oo 

g 2 (\k\/A coul ) = ( g 2 + J29 2 n+2 < )\ 90 = 90 m/* Q c D y ( 8 - 2 ) 

n=l 

so the new coupling constant agrees with go in lowest order. Thus it is a regular redefinition 
of the coupling constant, and g may be used for perturbative expansions. 

Corresponding to the invariant charge is an invariant ,9-function defined by 

m = \k\ dgm, ^ oul) - (8.3) 

d\k\ 



It may be calculated from 



P(g) = ^\k\ dg ° m/ ^ QCD) 

990 d \ k \ (8.4) 

dg 

P(9) = Po(9o)\ 9o (g), 
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where g (g) is obtained by inverting Eq. (8.2). This may be simplified by using 

- — — ' si.a Mm) (8.5) 



|fc|=A' 



2gx dg 

P AdV 
2gx dA 

by Eq. (6.2). From Eq. (4.8), this gives 

1 oo 

M = - Tq ( 2^ 4 + E<i^o n+2 )U), (8-6) 

where v' nl = v n i/xo, and we have used vu/xo = 26 - To find V(|fc|) to arbitrary accuracy, 
one calculates (3(g) perturbatively and then solves the flow equation (8.3) for g(\k\/A cou \). 

Finally we remark that we may choose new unrenormalized and renormalized expan- 
sion parameters according to 

9o = #(A/A CO ui) / 

(8-7) 

g' r = g(p/A cou i). 

so g' and g' r lie on the same invariant trajectory, the only difference between them being 
the value of the argument. 



9. Conclusion 

We have successfully applied the perturbative renormalization group to the Coulomb 
gauge. In this gauge, the 44-component of the gluon propagator D44(|fc|, &4, Aqcd) is a 
renormalization-group invariant in the sense that it is independent of the regularization 
and renormalization schemes, and of the ultra-violet cut-off A and renormalization mass 
\i. With the help of the perturbative renormalization-group we have decomposed it into 
an instantaneous part V(|fc|), which we call the color-Coulomb potential, and a vacuum 
polarization part P(\k\, k^) which vanishes at large k±. Each of these terms is separately a 
renormalization-group invariant, and their asymptotic form, at large \k\ and k<± respectively, 
was reported in the Introduction. 

The color-Coulomb potential allows us to define an invariant QCD charge g(\k\/ A cou \) 
by k 2 V(k) = xo(7 2 (|A;|/A cou i), where xo = n^ 2 ^ Nf ■ This invariant charge is the QCD 
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analog of the invariant charge of Gell-Mann and Low in QED. We have shown how to 
calculate the corresponding invariant /3-function, = /3(g). Because this charge is 

o\k\ 

scheme-independent it may offer some advantage in providing a definition of a s (k 2 ) = 
g (|fc|/A cou i) ^ w ] iereas t he standard definition in current use is scheme-dependent (see for 
example [11]). The color-Coulomb potential V(|fc|) is also the natural starting point for 
calculations of bound-states such as heavy quarkonium. 

The Coulomb gauge provides direct access to quantities of non-perturbative inter- 
est. Indeed both V(|fc|) and -P(|fc|, /C4) have a natural role in a confinement scenario: 
V(\k\, Aqcd) is long-range, anti-screening, and responsible for confinement of color-charge, 
whereas the vacuum polarization term P(\k\, k^) is screening, and responsible for "breaking 
of the string" between external quarks, when dynamical quark pairs are produced from 
the vacuum. We expect the linear rise (or not) of V(|if|) at large \x\ to provide an order 
parameter for confinement of color charge, even in the presence of dynamical quarks when 
the Wilson loop cannot serve this purpose. The accompanying article [1] reports a nu- 
merical study of the running coupling constant g 2 (k / A. cou \) . The data show a significant 
enhancement at low \k\, in agreement with this confinement scenario. However additional 
studies at larger values of (3 = 4/<7q are necessary before a conclusion can be reached about 
a linear rise of V(x) at large \x\ in the continuum limit, (3 — > 00. 

The data also show a strong suppression of the equal-time, 3-dimensionally transverse, 
would-be physical, gluon propagator D\ r Ak) at k = 0, and agree with a formula of Gribov 
that vanishes like \k\ near \k\ = 0. The only explanation for this counter-intuitive behavior 
is the suppression of configurations outside the Gribov horizon in the minimal Coulomb 
gauge. Since D\ v Ak) is strongly suppressed, we may understand the main long-range forces 
between color charge as being due to D44 which, as we have seen, is the sum of the attrac- 
tive instantaneous color-Coulomb potential V(|fc|) that is anti-screening, and the vacuum 
polarization term P(\k\ 7 £4) that is screening. According to the confinement scenario dis- 
cussed in Sect. 3 and in [2] and [1], both are long-range in the minimal Coulomb gauge 
because entropy favors a high density population close to the Gribov horizon. 
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Appendix A. Relation of Faddeev-Popov and phase-space functional integrals 

We wish to derive the canonical or phase-space functional integral (2.3) from the 
Faddeev-Popov formula (2.1). The argument merely reverses the text-book derivation of 
the Faddeev-Popov formula from the canonical phase-space functional integral in Coulomb 
gauge while keeping track of the sources J^. We introduce the identity 

exp[-(l/2) J d 4 xF^ }=N Jd 3 E exp J d 4 x[iE l F 0l - (1/2) E?], (A.l) 

which is a Gaussian integral over new variables Ef that will play the role of independent 
color-electric field variables. This allows us to rewrite (2.1) as 



Z(J, J 4 ) = / drA d E S(diAi) det(-A<%) 

Jg (A.2) 
x exp / d 4 x [iEi(Ai - D,A 4 ) - (1/2) (Ef + Bf) - igoJ^)), 

where Bf = — d^A^ + g§f ahc A\A%, etc. Integration on A 4 imposes color-Gauss's law, 
DiEi = go J a, in the form of the constraint S(DiEi — goJ^), 



Z(J, J 4 ) = / d 6 A d 6 E 6(diAi) det(-A^) 5(D t E t - g J±) 

(A.3) 

x exp J d 4 x [iEiAf - (l/2)(Ef + Bf) - ig JiA?)}. 

The constraint expressed by 5(diAi) has allowed us to replace Ai by its transverse part 
A\ T everywhere. We separate the transverse and longitudinal parts of Ei = Ef — di(f>, and 
we have d s E = NdE tr d(p. The Faddeev-Popov determinant is absorbed by 

det(-A^) 8(DiEi - g J 4 ) = det(M) 8(M(j) - p coul - g J±) 

(A.4) 

= S[(j)-M Vcoul+SO^)], 

where the symbols are defined as in Eq. (2.7). We now integrate over 0, and in a similar 
way we integrate out the longitudinal part of Ai, to obtain (2.3). 



Appendix B. One-loop expansion 

In this Appendix we find the one-loop expansion of the quantities Vq and Po defined 
in Eqs. (2.14) and (2.15), and which appear in -Do, 44 = Vq + Po. 
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The Faddeev-Popov operator is written M(A tr ) = M + M l (A tr ), where M = -df is 
the negative of the Laplacian, and (Mi) ac = —g f abc A b,tr di. The color-Coulomb potential 
energy functional V(A tr ), defined in Eq. (2.11), reads 

V(x,y;A tT ) = [(M + M 1 )- 1 M (M + Mi)" 1 ]^, (B.l) 

and has the expansion 

V(x, y; A tr ) = [M~ l - 2M~ 1 M 1 M~ 1 + 3M~ 1 M l M~ 1 M 1 M- 1 + . . ]^ , (B.2) 

where M~ x \ s ^ = (2tt)- 3 / d 3 k (k 2 )' 1 exp[ik ■ (x - y)} = (4tv\x - From Eq. (2.14) 

for Vq we obtain to one-loop order 

V (x -y) = gliM- 1 + 3M^( M l M~ 1 M l ) M" 1 ]^ S(x 4 - y 4 ), (B.3) 

where we have used (M 1 (A tr )) = 0, which holds because Mi(A tr ) is linear in A tr . The 
average designated by (. . .)o refers to the free-field average, with free-field propagators 
given in Eq. (3.1). This gives 

V = ^o,o + V ,i (B.4) 
where the zero-loop piece is given explicitly by 

V ,o(x ~ y)$ ae = 9 2 M^(x-y) 5(x 4 - y 4 ) 5 ae (B.5) 
and the one-loop piece by 

V 0>1 (x ~ y)* ae = %o / d 3 x'd 3 y' f abc f cde (A?> b (x',x 4 )Aj> d (y',x 4 ))o 

M^(x-x') d.M^ix' - y 1 ) djM^iy 1 - y) 5(x 4 - y 4 ). 

(B.6) 

These terms are illustrated in Fig. la. In momentum space we have Vb,o = 9o/k 2 : and 

V rifkl) = ^V)- 4 f d* P k * S «-KM k J . (B.7) 

' Ul ^ (k 2 ) 2 " ' J (p2 +p 2 A){p _ k) 2 ^ > 

The result of this integral is given in (4.2) and (4.3). 

Similarly, for P given in Eq. (2.15), we have to one-loop order 

P (x - y)8 ad = -g 2 ( (M -V c a ou i)Or) (M^p d coul )(y) ) . (B.8) 
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where p a coul = -g f abc Af' b E^ c . This gives 

P (x - y)5 ad = -sSf d3x 'd 3 y' M o _1 0? - *) 

x (p a coul (x\x 4 )p d coul (y',y 4 ))o M^tf-y), 

where 

(Pcoui(^)Pcoui(?/))o = 9 2 of abc f de9 [(Af\x)A^ e (y)) (Ef > c (x)Ef>° (y)) 

HAf\x)E^{y)) Q {E^{x)Af e {y)) Q ] 

This term is illustrated in Fig. lb. In momentum space it is given by 



(B.9) 



(B.10) 



-9oN A f 4 Pijip) Pijip-k) 2 



' {k 2 y y J " (p 2 + pD [(jr_ fc)2 + (p 4 - fc 4 ) 2 r 

(B.11) 

where (p) = 5^ — is the 3-dimensionally transverse projector. The contraction in 
the numerator gives 2 terms, 



P %J {p)P %J {p-k) = Jt + J 2 

(B.12) 



Ji = 2 



£2 fc 2 _ fc )2 
2 = _ w w ~ ' 

p 2 (p — k) 2 

Each term results in a Feynman integral I\ and /2- The integral I2 looks more complicated. 
However it is only logarithmically divergent by power counting, and when the integration 
is performed, the coefficient of the divergent part of I2 vanishes, so I2 is finite. As a result 
J 2 (|^|, & 4 ) vanishes in the limit fc 4 — > 00, and does not contribute to Pq s (/c). The result of 
the 1\ integration is given in (4.2) and (4.4). 

The integrals (B.7) and (B.ll) are evaluated by dimensional regularization, withp 4 — > 
Pd, and p = (pi) for % = 1, . . . , (d - 1). 
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Figure Captions 

Fig. 1. Horizontal lines correspond to instantaneous propagators 1/k 2 . Curved lines 
correspond to non-instantaneous propagators l/(k 2 + k 2 ). Diagram la is a con- 
tribution to Vq. Diagram lb is a contribution to Pq. 
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(b) 



Figure 1 
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